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Abstract 

It is shown that in realistic SUSY GUT models of quark and lep- 
ton masses both the proton decay rate and branching ratios differ 
in general from those predicted in the minimal SU(5) supersymmetric 
model. The observation of proton decay, and in particular the branch- 
ing ratio B[(p — > ir + v)/(p — > K + v)], would thus allow decisive tests 
of these fermion mass schemes. It is shown that the charged lepton 
decay modes p — ► K°{i + ,p — ► K°e + arising through gluino dressing 
diagrams are significant and comparable to the neutrino modes in 
large tan/3 models. Moreover, it is found that in certain classes of 
models the Higgsino-mediated proton decay amplitudes are propor- 
tional to a model-dependent group-theoretical factor which in some 
cases can be quite small. There is thus a natural suppression mech- 
anism which can explain without adjustment of parameters why in 
the context of SUSY GUTs proton decay has not yet been seen. The 
most interesting such class consists of SO (10) models in which the 
dominant flavor-symmetric contribution to the up-quark mass matrix 
comes from an effective operator of the form 16^16^10^45 jj, where 
(45//) points approximately in the I^r direction. This class includes 
a recent model of quark and lepton masses proposed by the authors. 
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1 Introduction 



It is well-known that in supersymmetric grand unified theories (SUSY GUTs) 
the dominant contribution to proton decay is through dimension-five opera- 
tors which arise from the exchange of superheavy, colored Higgsinos. 1,2 ' 3 In 
this paper we shall take a closer look at these decays. It will be shown that 
two simplifications that are usually introduced in the analysis of this effect 
are in certain interesting cases not justified, and that in going beyond these 
simplifications some important features emerge. 

The first simplification usually made is to analyze Higgsino-mediated pro- 
ton decay in the context of the minimal SUSY SU (5) model. In that model 
there is a very simple relation between the proton-decay amplitudes and the 
mass matrices of the quarks and leptons that allows one, in fact, to write these 
amplitudes in terms of the quark masses and the Kobayashi-Maskawa matrix 
elements. This relation comes from the fact that the Yukawa couplings of the 
colored Higgsinos which mediate proton decay are equal, because of SU(5) 
invariance, to the Yukawa couplings of the ordinary light Higgs doublets, H 
and H', that give rise to light fermion masses. But this cannot be taken 
seriously precisely because the minimal SUSY SU(5) model gives a wrong 
account of the light fermion masses. In particular, minimal SU(5) predicts 
that m" = m° d and m° = m° s (superscript zeros refer to quantities evaluated 
at the GUT scale), and thus m e /m^ = md/m s , which is off by an order of 
magnitude. 

The second simplification is to neglect gluino loops in dressing the ef- 
fective, dimension-five, AB ^ operators to make four-fermion operators. 
This implies that only W-ino loops need be considered, which in turn leads to 
the conclusion that the only significant modes are those that involve neutral 
leptons (p — > K + V and p — > -n + u). It will be explained why for models with 
large tan/5 the gluino loop diagrams cannot be neglected and decay modes 
with charged leptons (p — > K°fi + , p — > K°e + ) can become comparable to 
the neutrino modes. 4 (Since the Higgsino-mediated proton-decay rate goes 
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as tan 2 (3 it might be thought that the case of large tan (3 is in conflict with 
present limits. While this may be true in minimal SUSY SU (5) model, there 
are several likely suppression mechanisms which have been proposed in the 
literature 5 in going beyond minimal SU(5) and one such group theoretical 
suppression mechanism is found in section 4 of this paper. Large tan (3 is ac- 
tually a natural feature of many grand unified models, especially those based 
on SO (10) or larger groups where it often is predicted to be m t /m b .) 

Realistic models of quark and lepton masses in the context of SUSY GUTs 
require that the "bad" SU(5) or SO (10) mass relations be broken. This 
means that the quark and lepton masses receive contributions from additional 
operators not present in the minimal SU(5) model, operators which either 
involve Higgs fields in representations larger than the fundamental or are of 
higher dimension. What this implies is that the proton-decay amplitudes are 
related in a less direct, model-dependent way to the quark and lepton mass 
matrices. This allows one in principle to distinguish different theoretical 
models of quark and lepton masses by their proton-decay branching ratios, 
as will be seen. It will also be seen that an interesting group-theoretical 
mechanism exists by which the higgsino- mediated proton decay rate may 
be suppressed to levels near but consistent with present limits without any 
special adjustment of parameters. Interestingly, the group structure required 
for this suppression is precisely that suggested in a recent model of quark and 
lepton masses. 6 

2 Review of p Decay in Minimal SU(5) 
(a) W-ino loops: 

In order to set the stage for the later analysis it is convenient to review 
briefly the standard analysis of Higgsino-mediated proton decay in the mini- 
mal SUSY SU (5) model. In that model the Yukawa couplings of the quarks 
and leptons come from the following terms in the superpotential. 
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^Yukawa = \ [10,10;] 5 H /v + £ D tJ 

i,j=l i,j=l 



5,10, 



where C/jj and are the mass matrices of the charge- 1 and charge- 1 quarks. 
In terms of multiplets of the standard model gauge group this gives 



^Yukawa(^, = E ^ \<Q 3 \ H / V + E [<%Qj + Ulf\ H'/v', (2) 

ij 

and 

W / Yukawa(-ffc, #C?) = E 



Hc/v+^Dij hQj + dXj H' c /v 



ij 



(3) 

From the diagrams of Fig. 1 there arise two types of B-violating quartic 
terms involving the interactions in eq. 3. These are given by 



Mcvv 



ki 



ij kl 



(4) 



where Mq is the mass of the superheavy color-triplet Higgsino. Here the 
contractions of 377(2) L and SU (3) c indices are as follows: 



{LiQ^iQkQi) = e a ^(ud* - e^)(4d] - d p k u]) 
(dX 3 )Klt) = e aM {d?uf){v%lt) . 

In terms of component fields the operators in eq. 4 are of dimension five 
and contain two fermions and two scalars (squarks and/or sleptons). Since 
squarks and sleptons are heavier than the proton they must be converted 
into quarks and leptons by a gluino, W-ino or neutralino. This means that 
the diagram in Fig. 1 must be dressed by a vertex loop. 

It is generally argued that the gluino loops can be neglected. This argu- 
ment can be stated in the following way. The products of superfields QjQkQi 
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and d^u C jU% appearing in eq. 4 are totally antisymmetric in color. Since they 
are bosonic, they must also be totally antisymmetric in flavor. (Note: flavor, 
not family.) Now suppose that one neglects all flavor-dependence introduced 
by the gluino loop. Then, in the resulting four-fermion operators, the prod- 
ucts of quark fields, ipQ^Qk^Qi and V^V^V^j will have the original flavor 
structure of eq. 4, namely they will be totally antisymmetric in flavor. But 
since these are now products of fermions their antisymmetry in color and 
flavor implies that they must also be totally antisymmetric in spin, which is 
clearly impossible for three spin-| objects. 

The gluino loops, therefore, are only important to the extent that they 
introduce some additional flavor dependence. This can happen through the 
nondegeneracy of the different flavors of squark. But that gives a GIM-like 
suppression factor which is known to be very small, especially for the squarks 
of the first two families, because of the limits on squark non-degeneracy from 
the neutral Kaon system. Mixing in the squark/slepton mass-squared matrix 
is another source of flavor dependence. (In addition, there is a suppression 
factor of m u that comes from the fact that the four-fermion operator must 
involve only u rather than c or t quarks. But this is not an extra suppression, 
because including flavor mixing in the squark masses, this m u factor will be 
replaced by m c or m t .) What shall be seen in section 2(b), however, is that 
if tan /3 is large flavor-change arises in the squark mass matrix which, though 
very small, is nevertheless enough to make the gluino-loop contributions to 
proton decay significant and comparable to the W-ino contribution. 

Proceeding with the standard analysis, one concludes that the dominant 
contribution to proton decay comes from dressing the effective, dimension-five 
operators coming from eq. 4 with W-ino loops. (The photino and zino loops 
suffer the same problems as the gluino loops, with an additional suppression 
of a/a s . The Higgsino loops have suppressions of small Yukawa couplings.) 
Because the second operator in the brackets in eq. 4 contains only SU(2)l- 
singlet fields it can therefore be neglected. 

From the form of the first term in eq. 4 it is readily seen that the W-ino 
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loop leads to the following effective four-fermion operators: 



f 1 U2 T~) TJ c 

~ M c vv' 4n ly j iUklta l 3 'y 

x[«rff)(^.)(/(4,^) + /(4,^)) 

+(dtuf)(u]l-)(f(u k ,d l ) + f(d i ,u j )) (5) 
+(d% U P)(d]v j )(f(u k ,d i ) + f(u l ,l-)) 

+(u- k df)(u]l-)(f(d k ,u i ) + f(d l ,u j ))\. 

Here the color indices (a, (3, 7) and flavor indices (i, j, k, I) are understood 
to be summed over, and the fermion fields paired together in parentheses 
are spin-contracted to singlets. / is a loop integral with dimension of M _1 
defined by f(a, b) = ^ (_^l n ^ - [ a _ b ]\ For m, ~ m~ b » 

a b \ a W W J 

mw, f(a,b) ~ and for m~ a ~ mg < m^, /(a, 6) ~ ^-ln^. While 
deriving eq. (5), it has been assumed that flavor mixing in the squark sector 
is negligible, which should hold to a good approximation in supergravity 
models. We have verified that the effect of flavor mixing that occurs in large 
tan/5 models (see discussion in section 2.(b)) is negligible in so far as the 
W-ino dressing graphs are concerned. 

If one neglects contributions proportional to m u , then it is easily seen 
that only the third term in the brackets in eq. 5 contributes to proton decay. 
One can rewrite that term in the physical basis of the charge- 1 quarks, where 
Uij is diagonal, as follows 

3 3 

^«EE DflUuteafrffiviXdZvi), (6) 

j=l k=2 

where we have taken the / factors to be flavor independent (which should be 
a very good approximation since the squarks (and sleptons) must be nearly 
degenerate.) Here d\ = J2j{^KM)ijdj- The important conclusion is immediate 
that only the neutrino modes of proton decay are important. 
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Note the (paradoxical) result that the ratio of d to s quarks, which deter- 
mines the branching ratio of 7r to K modes, is controlled by the elements of 
the matrix Uij. In particular, the ratio of the coefficients of the (du)(du) and 
[(su)(dv) + (du)(sv)] terms is (neglecting terms proportional to m u ) given by 



C d d ^ m c sin 2 9 C + A s m t V t 2 d 



C ds -m c sin 9 C - A s m t VuV cb 



Or: 



where a = T/ Vtd a , and b = Vti > mtAs are quantities that can be measured in 
low-energy experiments. (Except for the phase of b. m c /m t has in general 
a non-trivial phase which cannot be measured in low energy experiments, 
but which can be predicted in specific models of quark and lepton masses.) 
A s is a short distance renormalization factor, proportional to the Yukawa 
couplings. If tan/3 is small (tan/3 < 10), A s — (1 — Y t /Yf) 1 ^ 12 where Y t = h\ 
is the top-quark Yukawa coupling-squared and Yf is its 'true' fixed point 
values, i.e., the value of Y t at the weak scale if Y t were infinite at the GUT 
scale. (Yf ~ 1.29 corresponding to a s (Mz) = 0.12). For tan/3 = m t /mb, 
A s — 1. If tan/3 takes intermediate values, A s is to be evaluated numerically. 
The parameters a and b are of order unity, present experiments allow a factor 
of 2 uncertainty in their numerical values. 



(b) Gluino loops, large tan/3 and charged lepton modes: 

The coefficient of the effective dimension-five operator, LiQjQkQi, in eq. 4 
depends on the mass matrix of the charge- 1 quarks, Uij. Thus, a gluino loop 
which dressed this operator would be correspondingly enhanced by containing 
a virtual t squark. Such a diagram is shown in Fig. 2. Since, however, the 
external charge-| quarks must be w's, for this to happen a flavor change must 
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occur in the gluino loop, such as is indicated by the t*u squark mass insertion 
shown in Fig. 2. This gives a compensating suppression; the question is by 
how much. If large enough flavor change can occur in the squark mass matrix, 
then the gluino loops can be important. As will now be explained, this is 
indeed the case for large tan ft. 4 

If one neglects the Yukawa couplings of the lighter two generations then 
there are only two Yukawa terms left. These can be written 



--Yukawa 




(9) 



Here and throughout, the unprimed quark fields denote the mass eigenstates, 
while t' denotes the SU(2)l partner of bi, and b' denotes the SU(2)l partner 
oft L . Thus *' = V* b t + V; b c + V* b u and b' = V tb b + V ts s + V td d. 

Corresponding to eq. 9 there are soft SUSY-breaking trilinear terms 



^-&{\H£)+*HIH%Y (l0) 

Here and throughout, the squark fields t, b, c etc. are not the squark-mass 
eigenstates but the superpartners of the quark-mass eigenstates t, b, c etc. 
Similarly, t' and b' are the superpartners of t' and b'. 

The interactions given in eq. 10 give rise to the diagrams shown in Fig. 
3. There are also related diagrams in which the vertices are the hard, SUSY- 
invariant Yukawa couplings and the SUSY breaking comes from the non- 
equality of the squark and quark masses. Altogether, these give effective 
low-energy terms of the form 



AC ^m 2 (3 + A 2 )^ 2 )' 1 ^ 



M. 



GUT 



SUSY 



Cl 



v' J 



rt' + c 2 



^)b>*b> 

V J 



(11) 
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where c\ and C2 are numbers of order one that depend on the precise spec- 
trum of sparticles. eq. (11) with c\ = C2 = 1 is what one would obtain 
with supergravity boundary conditions at the GUT scale for the squark fla- 
vor mixing, if one ignores the running of all the relevant parameters. This 
is, of course, very naive, but this should still tell us the correct results to 
within a factor of 3 or so. C\ and C2 parametrizes the effect of such running. 
Substituting t' = V t * b t + V* b c+V* b u, one sees that, for example, a u-t mixing of 
order (j^j Kb is produced. For large tan/?, ~ 1 and such flavor mixing 
becomes significant. 

Armed with this result one can go back and evaluate the graph shown in 
Fig. 2. One obtains (in the minimal SU (5) model) 



cfro „ 1 _^_ E DaU33KbVtj [ (zru _ Uif t ){djU) ] . (12) 

This is to be compared to the dominant W-ino graph contribution which 
gives 

C»~^EAii' I! [(»'«K»V,)]. (13) 

Since U33/U22 = ntf/m® one sees that for m^jv' ~ 1 (that is, for tan/5 ~ 
m>t/ m b) the gluino and W-ino contributions are comparable, and therefore 
the charged lepton and neutrino modes are comparable. 

To be more precise, we shall present below the effective Lagrangian aris- 
ing through gluino dressing to lowest order in the flavor mixing parameters 
= KU)y/[("i«,d)« - K d y for i ^ j for the case of minimal SUSY 
SU(5): 

+ «df (tijc) Aj (f(u h d>) - f( Ui , 4)) 
+ (df? «f)(«7e,)A]^ (f(u k , d' t ) - f( Uj , <)) 
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/(«*,<$■)) 

f(Uj,Uk)) 

fM)) 

fiUi.d'j)) 

M,<$)) 
f(uj,d' k )) 

/(«*<?■))]■ (14) 

In the above, we have used a basis where the up-quark mass matrix is di- 
agonal and the down-quark mass matrix is given by VkmD, D being the 
diagonal down mass matrix. It should be noted that in the approximation of 
keeping only the third family Yukawa couplings, the mixing parameter Af ■ 
in this basis is nonvanishing only for i = j = 3 because the only mixing is 
the down squark sector is of the form b'*b'. It is clear from eq. (14) that 
the dominant contribution for proton decay via gluino exchange arises from 
the first and the fifth terms, all other terms being proportional to the small 
up-quark mass. eq. (14) justifies the qualitative discussion preceeding it. 
In particular, if the SUSY spectrum is such that the gluino is lighter than 
the squarks, since rrig ~ (ctg/a^ra^, the glu* 110 dressing effect has a fac- 
tor of 4/3(a s /a2) 2 ~ 18 in the amplitude relative to W-ino dressing. Using 
the naive estimate of Ci = 1 in eq. (11) and A = 1, we find that the two 
amplitudes are about the same. 

It should be noted that, just as for the W-ino loops, the second operator 
in eq. 4, namely the d c u c u c l + operator, can be neglected in the gluino loops. 
This is because there is no flavor- changing in the right-handed squark masses 
coming from diagrams analogous to Fig. 3. If one interchanges t' and b c in 
Fig. 3a, for example, one merely gets the effective flavor-conserving operator 
b c *b c . 



+ (dfu f3 k )(u]e l )Af j (f(u k: d' i )- 

- «^)«^)A^.(/K,^)- 

- (uf<tf (/(«„<)- 

+ (rfrrff)(^)A^(/K,o- 

+ (^f)K^Aj (/(«) - 

+ «<)«^)Aj(/(^,4)- 
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3 Models With Realistic Fermion Masses 



So far we have been considering the unrealistic case where the quark and 
lepton Yukawa couplings are those of minimal SUSY SU(5). As we have 
seen, what controls proton decay in the general case (for both W-ino and 
gluino loops) is the Yukawa coupling of the colored Higgsino, H c , to the 
quark bilinear QiQj and of the H' c to LiQj. For minimal SUSY SU(5) 
the former is given by the up quark mass matrix, Uij, as shown in eq. 3, 
but in general it will be some different matrix, which we will denote U[j. 
W / Yukawa(^c) = Eij U-j [QiQj]H c /v + .... Similarly, the coupling of H' c is 
given by some matrix D'^ which is not in general the same as Dij. Let us 
now examine the relationship between U[j and U^ in some simple cases. 

The matrix U^ can be written as the sum of a symmetric and an anti- 
symmetric piece: 



Let us assume for simplicity that there is only a single term that contributes 
to the symmetric piece. Then the matrix U[, will have the form 



It is easily seen that the antisymmetric piece does not contribute to U[j, 
for € al 3jQfQjHQ being antisymmetric under SU{2) L and SU(3) C has to be 
symmetric in flavor. The factor r is group-theoretical and gives the ratio of 
the coupling of the color-triplet Higgs(ino) to the doublet Higgs(ino). This 
factor will be of great interest to us later. 

To see the branching ratios in proton decay we want to go to the physical 
basis of the charge- 1 quarks, as in the discussion leading to eq. 6. This is 
done by some transformation under which Qi — > WijQj, u\ — > WijUp and 
U — > W T UW, so that U' gets transformed as U' — > W T U'W. It is interesting 
to see what happens if there is a hierarchy among the elements of U, as is 



Uij = f(tj) +g[ij]. 



(15) 



U'ij = r- f(ij). 



(16) 
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usually the case in models of quark and lepton masses. Assume, therefore 
that > f'23 ~ #23 > f'22, and that the elements of the first row and column 
can be neglected. Then one has, approximately, that 



U 



( 





22 







/33 



while 



U' 



( 






22 









/23~^23 
/33 



\ 



/33 / 



\ 



( 







m c 




\ 


m t J 



(17) 




(18) 



923 

923 /33 / 

One sees that the relative contributions to proton decay from the (2,2) and 
(3,3) elements are the same as in the minimal SUSY SU(5) model. But 
now there is also the contribution from the (2,3) element, which is obviously 
model-dependent. In particular models of quark and lepton masses (723 would 
in general be known. If tan/3 <C m t /mi, then the W-ino loops dominate and 
proton decay is controlled by the terms in eq. 6 with the matrices D and U 
replaced by D 1 and U'. One can then generalize eqs. 7 and 8 using eq. 18. 



a 



dd 



— sin 9 r 



1 + a 2 b — 2ab2g 23 



(19) 



Cds \ 1 — ab — (a — l)b-2g 23 J 

where g 2 3 = 923/{ m c m <t)" 1 has been normalized so that it is typically a number 
of order unity (if non-zero). Thus 

T(p — > 7iV)/T(p — > Ku) 



[T(p -> 7vu)/T(p -> Ku)] min SU (5) 



1 - 2 



l +ab 2 I #23 



1 (o-l)W - 
1 i_„6 #23 



(20) 



One can straightforwardly derive the analogous expression when the first row 
and column of U' cannot be neglected. 
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What one sees is that by measuring the ir + V to K + V ratio in proton 
decay one could distinguish different models of quark and lepton masses, 
since in different models g 2 3 would take different — and for predictive models, 
computable — values. 

If tan (3 is large, then these neutrino modes can get comparable contribu- 
tions from both the W-ino and gluino loops. This would make the expression 
for the branching ratio T(p — > -K + u)/Y{p — > K + v) much more complicated 
(unless tan/5 is large and the spectrum is such that the gluino loop domi- 
nated). One can tell the relative importance of the W-ino loops and gluino 
loops by the relative importance of the neutrino and charged lepton modes, 
of course. 

A further test of models comes, if tan/3 is large, from examining the ratio 
of e + to This probes the value of D' n /D' 21 , as can be seen from eq. 12, 
where D and U should be replaced in the general case by D' and U' . In 
minimal SU(5), D'^ = D {j and T(p -> K°e + )/T(p -> fT /i + ) = \D n /D 21 \ 2 = 
| m d /m s sin 6 C \ 2 = 0.05. (Recall that we are working in the basis where Uij 
is diagonal, so that is nondiagonal.) In a realistic model of quark and 
lepton masses this branching ratio would be different. For example, in the 
Georgi-Jarlskog model 7 it comes out to be about 0.02. 

Finally, for large tan/5, the branching ratios T(p — > K°l + )/T(p — > tt°1 + ) 
give independent information about the elements of U[,-. (Cf. eq. 12.) 

4 Group-theoretical Suppression of p Decay 

Where eq. 16 applies one sees that the proton decay amplitudes coming from 
either W-ino or gluino loops are proportional to the group-theoretical factor 
r. This factor is given for various simple cases in Table I. 

In discussing the value of r it should be kept in mind that the form given 
in eq. 16, in which r appears, is only valid for cases in which the symmetric 
contribution to comes from a single term (or at least from several terms 
which all have the same gauge structure). The simplest possibility in SU(5) 
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is that the symmetric contribution to Uij comes from a 5 of Higgs, as in 
the minimal SU (5) model. In that case the couplings of doublet and triplet 
Higgs are trivially the same and r — 1. The other SU(5) multiplet that 
can couple flavor-symmetrically to the charge- 1 quarks is the 50 //. But, 
while this contains a color-triplet it does not contain the usual type of Higgs 
doublet, so that if only a 50 h were present then would vanish and r would 
be infinite, which is clearly unrealistic. If both a 5// and a 50// contribute 
then eq. 16 does not apply. 

One can also consider higher order operators involving bilinears of Higgs 
fields. Suppose that the flavor-symmetric contribution to comes from the 
operator [10^10^] 5 H 24 H . The product 5 x 24 contains 5 + 45 + 70. Of these, 
70 does not couple to 10 x 10, and 45 couples antisymmetrically. Thus there 
is only one term that couples symmetrically in flavor and eq. 16 does apply 
with a value r — — |, as can be straightforwardly shown. 

The next simplest possibility is that the symmetric contribution to 
comes from the bilinear 45// x 24//. There are two independent ways to 
contract this to get a flavor symmetric term, namely into a 5 or a 50, so eq. 
16 does not apply unless for some reason only a single contraction or linear 
combination of contractions of 45// x 24// contributes. For example, if only 
the operator [10jl0j] (45//24#) 5 appears, where 45// x 24// is contracted 
in the 5 channel, then r = 2/y / 3. (Such a contraction might come about 
by integrating out a 5 + 5, as shown in Fig. 4a. However, in that case 
there is also in general the direct 10jl0j5// term, which, if it contributes 
significantly to and U[p invalidates eq. 16.) A different term arises 
from integrating out a 10 + 10, as shown in Fig. 4b. If this gives the only 
symmetric contribution to then r = 0. (The reason for this is that the 
color triplet in the 45// does not have a coupling to a pair of quark doublets. 
See the discussion above of why does not contribute to U^.) This would 
be a way of suppressing Higgsino-mediated proton decay group-theoretically 
in SU(5). However, SU(5) seems to require an elaborate Higgs structure 
involving 75// + 50// + 50// to solve the doublet- triplet-splitting problem, 
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and the mechanism being discussed here would require in addition a 45 h- 
We thus consider the 50 (10) example of a group-theoretical suppression of 
proton decay discussed below to be more interesting. 

Turning now to 5*0(10), the simplest possibility is that the sole symmetric 
contribution to U^ is from a 10//, which, like the case of the fundamental 
Higgs in SU(5), gives r — 1. On the other hand, if the sole contribution to 
Uij is from a 126 h then r = y/3. (It should be noted that there are two 
color-triplets in the 126//, or more precisely two (3, 1, — |) representations 
of SU(3) C x SU(2) L x U(1) Y - [QiQj] couples to one linear combination of 
these, r = is the ratio of the strength of this coupling to the strength 
of the Higgs doublet's coupling to [u^Qj]. But the proton decay amplitude 
will also depend on the mixing angle which tells how much of the lightest 
color-triplet Higgsino is contained in the linear combination that couples to 

[QiQj]-) 

The simplest Higgs bilinear to be considered in 50 (10) is 10 H x 45//. 
This contains 10 + 120 + 320. Since the 320 does not couple to [16^16^-], 
and 120 couples antisymmetrically, there is only one symmetric contribu- 
tion from this term and eq. 16 applies. The value of r depends on what 
direction in group space (45//) points in. Let us call that direction, which 
is a linear combination of 50(10) generators, Q. Then it is easily seen that 
r = ^ q~+q ~ ) • (Note that the same expression applies to the SU (5) case 
of 5// x 24// with Q = Y/2 and gives r = 2(|)/(| - |) = -§.) This is made 
small if Q points approximately in the I$r direction, since the left-handed 
quark doublets are singlets under SU {2) R . One can write Q as a linear com- 
bination of I3R and Y/2. (There is a two-dimensional space of generators 
of 5O(10) that commute with the generators of SU(3) C x SU(2) L x U{l) Y .) 
Since Y/2 is non-zero for the left-handed quark doublets, it is obvious that 
r will be small only if Q points approximately in the I%r direction. 

This result is quite interesting, since, in a recent paper, 6 a model of quark 
and lepton masses is proposed by us in which the matrix U^ arises from an 
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effective term [16jl6j] 10#45// and the (45//) points approximately in the 
I 3 r direction. Indeed, in that model, Q ~ I 3 r explains three relations among 
known quantities: the Georgi-Jarlskog relation m° /m Q s ~ 3, the smallness 
of the second generation masses compared to those of the third generation, 
and the smallness of V c b- (In the "long version" of that model 4 Q ~ I 3 r 
also explains why m°/m° <C m°/mj.) What we have found here is that the 
same group theoretical assumption also gives rise to a (needed) suppression 
of Higgsino-mediated proton decay. 

In the model of Ref. 6 the group-theoretical suppression of the proton 
decay amplitude is numerically of order 10 _1 . This means that, while sup- 
pressed, p-decay is still within the reach of Super Kamiokande. 

It is worth noting that the term [16jl6j] 10//45// has another beautiful 
property, which is exploited in Ref. 6. For i = j = 3, 10 h x 45 # must 
be in the symmetric product, and, in particular, be contracted to into a 10, 
which gives equal contributions to D 33 and L 33 . But for i,j = 2, 3 or 3, 2 
contractions into both 10 and 120 are allowed, and the latter contributes 
differently to D and L. Thus a term of the form [16jl6j] 10 H 45 H explains 
why m° = m° while m° ^ m° . 

5 Conclusions 

We have shown that in realistic models of quark and lepton masses in the 
context of SUSY GUTs, the proton decay rate and branching ratios are 
different from the predictions of minimal SUSY 377(5) model. This would 
enable one to test quark and lepton mass schemes in SUSY GUTs by mea- 
suring, for example, the branching ratio B[(p — > 7i + V)/(p — > K + v)]. In 
predictive models of fermion masses, such branching ratios are computable 
as they are related to other low energy observables. We have also emphasized 
that in large tan/3 models, as often occurs in 5*0(10) GUT, the gluino dress- 
ing of the effective AB ^ dimension 5 operator is quite significant. This 
opens up the possibility that the charged lepton decay modes of the proton, 
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p — > K°fi + ,p — > K°e + , may be comparable to the neutrino modes. We give 
simple analytic arguments showing how this happens, especially in the case 
of tan/5 = m t /mi, schemes. The branching ratio B[(p — > K°fi + ) /{p — > fT°e + )] 
will then provide additional clue to the texture of the fermion mass matrices. 
The relative strengths of the neutral to charged lepton modes will tell us 
about the parameter tan/3 itself. While the precise values of these branch- 
ing ratios depend on poorly known (presently) Kobayashi-Maskawa mixing 
angles (V u b, V t d), that situation should change in the near future. 

We have also found, as discussed in section 4, an interesting group the- 
oretical suppression mechanism for proton decay rate. This could provide a 
simple answer to the question why in the context of SUSY GUTs, proton 
decay has not yet been observed. We find it interesting that a recent model 6 
of quark and lepton masses proposed on quite independent grounds automat- 
ically has this ingredient needed for group theoretical suppression of proton 
decay rate. 

In closing, let us emphasize that the discovery for proton decay would 
not only provide evidence for the violation of baryon number symmetry near 
the Planck scale, but it would also provide, in the context of SUSY GUTs, 
important clues to the structure of the fermion mass matrices. The search 
for these decays should continue in earnest. 
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Table I: The group-theoretical factor r that enters into the proton- 
decay amplitude in SUSY GUT models for which the flavor-symmetric con- 
tribution to the mass matrix of the up quarks comes from a single kind of 
operator. The second column gives the product of Higgs fields appearing 
in that operator. In the case of the 45# x 24h of SU(5) more than one 
contraction of indices is possible. 



Gauge group Higgs operator r 



SU(5) 


5h 


1 


SU(5) 


50 h 


oo 


SU{5) 


5 H x 24 H 


2 
3 


SU(5) 


(45 H x 24 H ) 5 


2/VS 


SU{5) 


(45 H x 24 H ) Fig . 2h 





50(10) 


10* 


1 


50(10) 


126^ 


V3 


50(10) 


10 h x 45 H 


Q(u,d)+Qu c 
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Figure Captions 



Fig. 1: Diagrams involving the exchange of supermassive, colored Higgsinos 
which give effective dimension-five AB ^ operators. The dominant 
contribution to p decay in general SUSY GUTs comes from (a). 

Fig. 2: A dimension-five AB ^ operator dressed by a gluino loop to 
produce a four-fermion interaction. By having a virtual t-squark this 
diagram is enhanced by a factor of U 33 = m t . This requires the flavor- 
changing squark-mass insertion indicated by the solid circle. 

Fig. 3: A one-loop diagram giving a flavor- changing squark mass term. 

Fig. 4: Two diagrams that would give operators of the form 10jl0j45#24#. 
Diagram (a) gives the 45# and 24# contracted into a 5. Diagram (b) 
gives a different contraction, for which r = 0. 
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Fig. la 




Fig. lb 
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Fig. 3b 
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